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■ ABSTRACT 

(N 

, Context. Intermediate mass galaxies (>10"' Mq) at z~0.6 are the likeliest progenitors of the present-day, numerous population of 

. spirals. There is growing evidence that they have evolved rapidly since the last 6 to 8 Gyr ago, and likely have formed a significant 

fraction of their stellar mass, often showing perturbed morphologies and kinematics. 
QQ ' Aims. We have gathered a representative sample of 88 such galaxies and have provided robust estimates of their gas phase metallicity. 

, Methods. For doing so, we have used moderate spectral resolution spectroscopy at VLT/F0RS2 with unprecedented high S/N al- 

lowing to remove biases coming from interstellar absorption lines and extinction to establish robust values of Rj-i = ([0//]/13727 -l- 
[0///]^^4959, 500T)/Hj3. 

" Results. We definitively confirm that the predominant population of z~ 0.6 starbursts and luminous IR galaxies (LIRGs) are on aver- 

r> ] age, two times less metal rich than the local galaxies at a given stellar mass. We do find that the metal abundance of the gaseous phase 

of galaxies is evolving linearly with time, from z=l to z=0 and after comparing with other studies, from z=3 to z=0. Combining our 
results with the reported evolution of the Tully Fisher relation, we do find that such an evolution requires that ~30% of the stellar 
mass of local galaxies have been formed through an external supply of gas, thus excluding the close box model. Distant starbursts & 
LIRGs have properties (metal abundance, star formation efficiency & morphologies) similar to those of local LIRGs. Their underlying 
physics is likely dominated by gas infall probably through merging or interactions. 

Conclusions. Our study further supports the rapid evolution of z~0.4 -1 galaxies. Gas exchanges between galaxies is likely the main 
cause of this evolution. 

Key words. Galaxies: evolution - Galaxies: ISM - Galaxies: spiral - Galaxies: starburst - Infrared: galaxies 



1. Introduction 

There is a growing consensus that most of the decline of star- 
formation density since z= 1 is related to the strong evolution in 
the intermed iate-mass galaxy population, with Mst ei from L5 to 
15xlOi°Mo (iHammer et alJl2005t iBell et alJl2005l) . At interme- 
diate redshift these galaxies are progenitors of present day spi- 



Send offprint requests to: myriam.rodrigues@obspm.fr 

* Intermediate MAss Galaxy Evolution Sequence, ESO programs 
174.B-0328(B), 174.B-0328(F), 174.B-0328(K), 66.A-0599(A) 



rals. This work belongs to a series of studies on intermediate- 
mass galaxies at z~0.6, in the frame of the large program 
IMAGE S (Interme diate MAss Galaxy Evolution Sequence). In 
Paper I (Yang et al.. 2008), we present the 2D kinematics of a 
representative sample of 65 galaxies and put into evidence the 
hig h fraction of non -relaxed galaxies at this redshift. In Paper II 
(Nei chel et al]l2008l) . we establish a morpho-kinematical classi- 
fication and find a strong evolution of the fraction of well-relaxed 
spiral rotating disk along the last 6 Gyrs. We also point out 
that intermediate-mass galaxies have double their stellar mass 
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since z~Q.6 according to the evolution of the TuUy-Fisher rela- 
tion found in Paper 3 dPuech et al.ll2008l) . As a whole, this se- 
ries of papers have shown that intermediate-mass galaxies have 
performed a strong evolution of their kinematical and morpho- 
logical properties since z~0.6. These results reveal the agitated 
history of the spiral disk galaxies over the last 6 Gyrs. 

The evolution of the metal content of the gas in galaxies is an 
useful tool to probe various galaxy evolution scenarios. In fact, 
the stellar mass-metallicity relation (M-Z) can help to disen- 
tangle the contribution of several processes taking place during 
galaxy evolution such as star-formation history, outflow pow- 
ered by supernovae or stellar winds and infall of gas by merger 
or secular accretion. In th e local Universe the M-Z relation has 
been w idely studied since iLequeux et aTl (ll979l) . lTremonti et alj 
( l2004 l) performed a reliable estimation of oxygen abundance and 
stellar mass for 53 000 star-forming galaxies from the SDSS. 
Several studies have characterized the relat i on at higher red- 
shift: < z < 1 (Kobulnic kv & Kewleyll2004t iLiang et^ai]|2006t 
Savagho et al. 2005; Lamar eille et al.ll2007h. l<z <2 dLiu etalJ 
2008; Erb et al. 2006; Shaplev et al.l2004l;lMaier et a l.'2006') and 
z>3 (Maiolino et al. 2008). A significant evolution is found as a 
function of redshift: z~0.7 galaxies show on average a gaseous 
abundance with a factor of two less than local counterparts at a 
given stellar mass jLiang et alj|2006h . Such a trend is also found 
at higher redshift. However, the metal abundance of the gas in 
galaxies is estimated for all these samples from strong line cali- 
bration as /?23 which uses [O iii], [O ii] and H/3 or N2, which uses 
[N ii] and Ha. Several parameters affect the measurement of line 
ratios: extinction and underlying Balmer absorption in the case 
of R23, the unknown contribution of primary nitrogen production 
as a function of the me tallicity for A^2, an d for both parameters 
the spectral resolution. ILiang et al.l (l2004b ) emphasized the ne- 
cessity of a spectral resolution over R > 1 000 and a S/N>10 to 
recover robust physical quantities from spectra. There are still 
few studies reaching the spectral quality necessary to estimate 
reliable metal abundance for intermediate redshift galaxies. The 
data from the Keck Redshift Survey are limited by the absence 
of flux calibration, making impossible to estimate the extinction. 
Moreover, metallicity estimated using the eq uivalent widths pro - 
duce systemati c ally higher metal abundance (Liang et"aLll2006l) . 
ISavaglio et al.l ( 12005,) observed 56 galaxies from the Gemini 
Deep Deep Survey with flux-calibrated spectra, but S/N of the 
data were insufficient to measure extinction for individual galax- 
ies. 

The aim of this paper is to obtain a robust M-Z relation for a 
reasonably large sample of intermediate mass galaxies at z~0.6. 
In this paper, we have followed the methodology proposed by 
ILiang etalJ (l2006) consisting of a precise estimate of the extinc- 
tion and underlying Balmer absorption in order to obtain reli- 
able /?23 abundance determination. The high quality spectra were 
taken with the VLT/FORS2. The paper is organized as follows: 
Sect. 2 describes the observational data and discusses the com- 
pleteness of the sample. In Sect. 3, we outline the methodology 
followed to estimate extinction, star-formation rates and metal- 
licities, as well as determinate the contribution of AGN. Results 
concerning the evolution of the M-Z relation and the effect of 
galaxy morphology are given in Sect. 4. In Sect. 5, we discuss 
systematic effects and the implication of our results on the close 
box scenario. We summarize our conclusions in Sect. 6. 

Throughout this paper we adopt a A-CDM cosmological 
model (Ho=70kms-iMpc-i, Qm=0.3 and Qa=0.7). All mag- 
nitudes used in this paper are in the AB system, unless ex- 
plicitly noted otherw ise. The adopted sola r abundance is 12 + 
log(0/H)^S.66 from lAsplund et all (|2004 . 
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Table 1. Number and type of objects in the redshift quality 
classes: secure (2) for spectra contain more than two strong fea- 
tures, insecure (1) when only one features was detected but with 
not very reliable supporting features, and single emission line (9) 
when the spectra had a single emission line without any features 



2. IMAGES observation 

2.1. VLT/F0RS2 Observations 

Galaxies were observed with the VLT/FORS2 in the Chandra 
Deep Field South (CDFS) in the context of the ESO large pro- 
gram "IMAGES". We randomly selected 270 galaxies in 4 fields 
of 6'8 X 6'8. Observations were performed in MXU mode, using 
two holographic grisms, GRIS 600RI+19 and GRIS 600z+23, at 
spectral resolution of 6.8 A and covering a maxima wavelength 
range from 5 120 A to 10 70 A . The slit masks w ere prepared 
using the GOODS catalogue ( iGiavahsco et al.l2004l) and FORS2 
pre-images. The arrangement of the slits has been optimized to 
improve the number of objects observed. 270 slits have been 
placed among which 241 objects have Iab <23.5 mag. The min- 
imum exposure time of each object in each grism is around 3 
hours. 

Spectra have been extracted and wavelength calibrated using the 
IRAFQ package. The spectra were flux calibrated using photo- 
metric standard stars. We have compared the spectrophotometry 
to HST/ACS photometry in V, I, Z bands and V, I, R in EIS cat- 
alogue and found a good agreement. 



2.2. Samples 

From the 241 target objects with Iab <23.5mag from the 
FORS2-IMAGES sample, we have identified 231 targets. For 
10 targets we have not been able to perform an identification 
because of instrumental problems (target falling outside the slit 
or bad sky extraction). Each redshift has been associated with 
a quality flag denoting the reli a bility of the identification, see 
description in iRavikumar et alJ (l2007h . We have classified ob- 
jects into 4 types: galaxy, split into emission line galaxies (ELG) 
and absorption line galaxies (ALG); quasar and stars. Table 1 
shows the number and type of objects in the different redshift 
quality classes. From the 270 slits of the FORS2-IMAGES, we 
have estimated the redshift for 240 objects, giving a success 
rate of ~88%. For objects with Iab <23.5mag the success rate 
reached ~96%. In FiglT] we have plotted the redshift histogram 
for all galaxies and QSO. The distribution have a range of [0.01- 
3.499] and a median value of 0.667. Two main peaks are visi- 
ble in the redshift distribution at z=0.670, z=0.735. These peaks 
correspond to s tructures in the CDFS field, already identified in 
IRavikumar etal.1 (120071) . 

From the Iab <23.5 mag objects we have selected 74 galax- 
ies having [On], H/?, [Oiii] emission lines, in order to evaluate 



' IRAF is distributed by the National Optical Astronomical 
Observatories, which is operated by the Association of Universities 
for Research in Astronomy,Inc., under cooperative agreement with the 
National Science Foundation. 
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Fig. 1. Redshift distribution (bin=0.02) of the FORS2-IMAGES 
sample, including the classes 2(secure), 1 (insecure) and 9(single 
line). 
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Table 2. Number of galaxies (208 galaxies + 4 QSO) by spec- 
tral type in the IMAGES and selected sample: [On], H/3, [Om] 
and H(j, emission lines detected (El); [On], H/3, [Oiii] detected 
(E2); [O ii] out of the wavelengh range (E3); Hfi or/and [O iii] out 
of range (E4); one or two emission lines detected (E5); galax- 
ies without emission lines (Abs). The E3 galaxies correspond to 
low-redshift objects. In this work we have only selected the 5 
galaxies with z>0.4 



their metallicities. Table 2 describe the construction of the work- 
ing sample by spectral lines availability. From the sample of 74 
galaxies, we have excluded 6 objects because of bad spectropho- 
tometry or wavelength calibration problems and 10 objects with 
H/3 or both [O iii] lines corrupted by strong sky emission or ab- 
sorption lines. Finally, our sample, hereafter sample A, is com- 
posed by 58 galaxies with a median redshift of 0.7. The S/N of 
H/3 line has a mean value of 40. 

The CDFS field has the advantage of being a widely studied 
field. The Great Observ atories Origins Deep Survey (GOODS) 
dGiavalisco et al.l |2004|) provides deep mu ltiwavelength obser- 
vations. HST/ACS and EIS photometries dArnouts et alj|2001h 
have been used to check the spectrophotometry and to calculate 
the aperture correction. The IR luminosity of galaxies has been 
estimated from the mid-I R catalogue of Le Fl oc'h et alj (l2005h 
and using the procedure of lCharv & E lbaz (20Ql]). Observations 
in X-ray from the Chandra X-Ray Observatory 1 Ms obser- 
vations (Giacconi et al. 2002) and radio detecti on from the 
Australia Telescope Compact Array 1.4 GHz dAfonso et alj 
I2OO61) have been used to identify Active Galactic Nucleus 



(AGN), see Sect 3.3. The rest frame magnitude in J,K-band 
and the luminosity at 2800A have been derived by modeling 
galaxy SED s using ISAAC and ACS multi-band photometry. 
iPuech et al] (|2008) have evaluated that the random uncertainty 
of this method is less than 0.2~mag in Mk- A catalogue of 
Sample A with RA, DEC, z, Iab from the ACS, M^, My, Lm 
is given in Table |6l 

We hav e comp leted our sample with 30 galaxies from 
iLiang et akl ( l2006h . referred hereafter as sample B. This 
Iab <22.5 mag galaxies have been randomly selected in Canada 
France Redshift Survey (CFRS) 3h, Ultra Deep Survey Rosat 
(UDSR) and Ultra Deep Survey FIRBACK (UDSF) and ob- 
served with VLT/FORS2/MXU. This sample ha s Mk measure- 
ment and mid-IR detection w ith ISOCAM, see iHammer et al.l 
(120051) and lLiang et al] d2004al) for a complete description. 

2.3. Completeness 

In this work, we have combined samples A and B and called 
the whole sample: sample A+B with a total of 88 galaxies. 
The two samples of galaxies span the same redshift range. In 
terms of absolute magnitude, sample B is composed of galax- 
ies brighter than sample A. However, the combination of this 
two samples gives a sample representative of the intermedi- 
ate mass galaxies at z~0.6. Fig. |2] shows the distribution of 
the K-band absolute magnitudes of sample A+B and the ob- 
jects with redshift identification in the F0RS2-IMAGES sam- 
ple verifying Iab <23.5mag. We have compared, on the same 
plot, the luminosity distri bution to the luminos ity function at 
redshift of 0.5 and I from iPozzetti et al] d2003h . Kolmogorov- 
Smirnov tKS tests support that sample A-i-B follows the LF at 
a 85% probability in the redshift range z=0.4 and z=0.98, and 
with MK< -21.5. The corresponding stellar mass range of com- 
pleteness is log(Ms,f//a,./MQ) > 10. The aim of the IMAGES 
large program is to investigate intermediate mass galaxies with 
M stellar > 1.5 X IO'^Mq, which represents 76% of our sample. 
Then, the incompleteness of the sample under Mk = -21.5 does 
not affect our conclusions. 



3. Methodology and data analysis 

The following section describes the methodology of data analy- 
si s for the sample A. The description for th e sample B is given 
in' Liang et all d2004al) and lLiang etal] d2006). The methodology 
used for sample A is very similar to the one used for sample B, 
thus the two sample are homogeneous. 

3. 1 . Flux measurement 

In order to derive metal abundance in galaxies with sufficient ac- 
curacy, it is necessary to estimate the extinction and the underly- 
ing Balmer absorption. We have corrected the Balmer emission 
lines from underlying stellar population. For each galaxy, we 
have fitted the observed spectrum, including its continuum and 
absorption lines, with a linear combination of ste llar hbraries. 
We ha ve used a set of 15 stellar spectra from the iJacobv et al.l 
(1 1 984 stellar library, including B to M stars (e.g B, A, F, 
G, K and M) with stellar metallicity. Spectra and stellar tem- 
plates have been degraded to the same spectral resolution of 8 A. 
The smoothing of the spectra emphasises the absorption fea- 
tures used in the fit. Only the continuum and absorption lines 
have been smoothed in the galaxy spectra. We have convolved 
the continuum, except at the location of emission lines, using 
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eratureof l OOOO K, the predicted ratio is 0.466 for I(Uy)/I(Ufi) 
lOsterbrockl 19891) . Using the following relation: 
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Fig. 2. Number counts (in logarithmic scale) of selected galaxies 
versus AB absolute magnitude in K-band. The black histogram 
refer s to sample A+B g alaxies. The luminosity function derived 
from iPozzetti et all (I2003.) at z = 0.50 is plotted in dashed line 
and at z = 1 in soUd line. The red vertical line represents the 
limit of 85%-completeness. 



the software developped by our group jHammer et al. I I2OOTI) . 
The best fit has been don e with the STARLIGH T software 
dCid Fernandes et al.ll2005l) . A ICardelli et all (1 19891 ) reddening 
law has been assumed. A rest-frame spectrum of one typical 
galaxy of our sample, J033210. 92-274722. 8, is given in Fig[3] 

We have chosen to use stellar templates rat her of the usually 
used SSP models of iBruzual & Chariot! ( |2003|) . hereafter BC03. 
Indeed, B C03 models overestimate the absorption in the H/3 line. 
lAsari et al . (2007); Cid Fernandes et al. (2005) have suggested 
that the origin of this biais is due to the STELIB library used in 
the BC03 models. The STELIB spectra have an excess of flux 
in both side of the H/3 when compare to model spectrum. As the 
synthetic spectrum is not use to retrieve physical properties of 
the stellar population, we have adopted a simple stellar library 
in order to better control the parameters of the fit. 

After subtracting the stellar component, we have measured 
the flux of emission lines using the SPLOT package. When the 
[Om] A4959 emission line was not detected, the flux was as- 
sumed to be 0.33 times the [Om]/i5007 (the ratio of the transi- 
tion probability). 

3.2. A reliable estimate of the extinction 

We have measured extinction by two methods in the sa mple A: 
Balmer decrement and IR/Hy8 energy balance (Lian g et"an 
l2004ah . Comparing the two values, we can check the reliability 
of our methodology and avoid systematic effects on abundances 
due to bias in extinction evaluation. 

For 37 objects of sample A, the Hy line has been detected 
and we could derive extinction in the gas phase directly from 
the spectra, using the observed Hy/HyS ratio. Assuming a case B 
recombination with a density of 100 cm"-' and an electronic tem- 



c(/(Hy)-/(H/3)) 



(1) 



we have determined the dust e xtinction c. Th e extinction param- 
eter Ay is calculated following|siato^ ([197^: Ay = E{B-V)R = 
(c^)/1.47 where R - 3.1. The median extinction of the sample 
is Ay - 1.53. 

Because of large uncertainties related to the measurement of 
the Hy line, we needed to verify the quality of our derived ex- 
tinction. Another method to evaluate Ay is from the infrared and 
optical SFR. Indeed, when dust extinction is ta ken into accoun t, 
the S FRu^ is in good agreement with S F/?ir (Flores et alj2004l) . 
The ratio of the uncorrected extinction S FRiia and S F/?ir allows 
to estimate the amount of extinction. The two SFR have been 
computed assuming the Kennicutt (1998) calibrations an d both 
use t he same Salpeter's initial mass function (IMF) (Salpete^ 
119551) . The Lir is estimated from th elCharv & Elbaz (2001.) re- 
lation and the 24 observations dLe Floc'h et al.ll2005h . The 
Liia has been estimated from the flux of Hj3 and assuming 
H„/H/3=2.87 (Osterbrock 1989). The Hfi flux have been cor- 
rected by an aperture factor derived from photometric magni- 
tudes at /ab and Vab and spectral magnitudes. Finally, the ratio 
has been corrected using the average interstellar extinction law: 



Ay(IR) = illog^:^ 
1.66 ^SFRu„, 



(2) 



The median extinction in the sample for lR/H/3 estimation is 
Ay = 1.71. In the case of the objects without IR detection, we 
have used the detection limit of 0.08 mJy of the MIPS catalogue 
to find an upper limit for Ljr and so an upper extinction limit 

Ay(IRli,n). 

In fig |4] we compare the extinction given by the two meth- 
ods. We plot in the same figure the objects from sample A+B 
with Ay (Balmer) and Ay(IR) measurements. The same method- 
ology has been used to determinate Ay(Balmer) and Ay(IR) in 
Sample B, see Liang et al. (2004a) for a complete description of 
the method. Most galaxies fall within the +0.64 rms 1-c r dis- 
persion. Th i s result is consistent with iFlores et"ail (12004 ) and 
iLiang et akl (l2004al) . We notice that the median extinction of 
the sample B is higher than sample A: Ay(Balmer)=1.82 and 
Ay(IR)=2.18. This difference between the two samples can be 
due to the selection of the sample B: galaxies are more luminous 
than in sample A and detected in near IR. 

Some discrepancies between the two measurements of ex- 
tinction are due to geometrical properties, as J033212.39- 
274353.6, J033232.13-275105.5 and CFRS 03.0932 which show 
Ay(IR) greater than Ay(Balmer). Indeed, these 3 objects are 
edge-on galaxies. In such case, a large fraction of the disk is hid- 
den by the dust in the disk plane. The detected optical Balmer 
lines just trace the star formation of a few optical-thin HII re- 
gions lying on the periphery of the galaxy in the line of sight of 
the observer. The consequence is an underestimated extinction 
value. The infrared radiation is less affected by dust so it comes 
from deeper region of the galaxy. Thus, we expect to have higher 
extinction values when estimated with the IR flux. 

For non-edge-on galaxies, we have assumed that the 
difference between the two measurements is statistical and 
not intrinsic: i.e. the two extinctions are statistically equal. 
Hereafter, we assume Ay as the mean value between Ay(IR) 
and Ay(Balmer) for the objects with IR detection and the mean 
between Ay(IRiiin) and Ay(Balmer) for galaxies with an IR 
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Fig. 3. Rest-frame spectrum of one of the sample A galaxies, J033210. 92-274722. 8. Grey boxes delimit the wavelength region 
where spectra may be affected by str ong sky hnes. Top panel: Spectrum with the location of strong emission lines. Bottom panel: 
Spectrum and synthetic spectra using lJacobv et al.l ([1984.) templates plotted in dashed line. The position of the absorption hnes are 
marked in dashed vertical hnes. 



upper limit. For objects without Hy detection we have used 
Ai/(IR) or Av(IRiim) when Ljr is not available. For edge-on 
galaxies, the difference between the two extinction estimates 
is physical as explained above. Since our aim is to correct the 
optical spectra, we have used Av(Balmer) for these galaxies. 
The comparison between the two extinction estimates allows to 
minimize systematic errors. In fact, Av(IR) tends to give higher 
extinction values because the mid-IR light is less affected by 
dust. On the other hand, the Balmer decrement may under- 
estimate the extinction because it follows only massive stars 
population and regions without dust-obscured star formation. 
Adopting the mean values between the two methods allows to 
reduce the systematic errors of each one. 



3.3. Searching for possible AGN contamination 

Before studying metallicities in galaxies we have to identify the 
objects whose lines are affected by contamination from an AGN. 
This identification is essential because the AGN processes affect 
the [O III] emission lines and so any metallicity estimate based in 
this line. First, we have eliminated from Sample A two galaxies 
harboring AGN spectral features, like broad Mg ii and Balmer 
hnes: J033208.66-274734.4 and J033230.22-274504.6. These 
two galaxies are very compact and are both X-ray detected. 
Broad Balmer lines suggest that these galaxies are Seyfert 1 
galaxies. 

Then, the diagnostic diagram log([Oii] /i3727/H^) vs 
log([Oiii] AA4959, 5007/H/j) was used to distinguish the HII 
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Fig. 4. Relation between the extinction Ay values derived by 
Balmer decrement and by SFRij^ and 5'F/?H„nc ratio. We plot 
the 37 galaxies from Sample A having reliable Hy detection and 
the 13 galaxies from Sample B with Ai/(Balmer) and Ay(IR) 
measurements. The objects with infrared detection are plotted in 
filled red circles. The arrows are the 21 objects in the Sample A 
with only Av(IR) upper limit estimation. The two dashed lines 
refer to the results with +0.64 rms. The down left image in the 
graph is the combining ACS/HST image in B,V and I band of 
the edge-on galaxy J033212.39-274353.6. 
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Fig. 5. Diagnostic diagram for the Sample A. Galaxies with 
only Ai/(IRiim) are plotted with down arrows. The objects with 
log([Neiii]^3868/[Oii]^3727) > -1.3 are plotted in red dots 
The so lid line shows the theoretical sequence from lMcCall et al.l 
(1 19851) for extra-galactic HII regions. The dashed Une shows 
the photo-ionization limit for a stellar temperature of 60 000 K 
and empirically delimits the Seyfert 2 and LINER from the H II 
regions. The dot-dashed li ne shows the dem arcation between 
Seyfert 2 and LINERs from lOsterbrockl (1 1 989t) . 
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region-like objects from LINERs and Seyfert galaxies, see 
Fig.|5] We have found 10 objects in the LINER and Seyfert 2 
region of the excitation diagram: J033214. 48-274320. 1, 
J033219.32-274514.0, J033222.13-274344.5, J033223.06- 
274226.3, J033224.60-274428.1, J033229.32-275 155.4, 
J033236.72-274406.4, J033240.04-274418.6, J033243.96- 
274503.5, J033245.5 1-27503 1.0. From these objects we discard 
5 galaxies which fall out of the H II region, even if we account 
for error bars: J033219.32-274514.0, J033222. 13-274344.5, 
J033240.04-274418.6, J033243.96-274503.5, J033245.51- 
275031.0. The other 5 galaxies are very close to the limit of 
star-forming galaxies to classify them definitely as AGN. We 
have chosen to keep them and we have used a different symbol 
for these objects in the figures. 

We search for evidence of shock processes in our galaxy 
sample by the presence of the emission line [Ne Iii]/13868. Five 
galaxies present log([Neiii]/13868/[Oii].13727) > -1.3 and fall in 
the Seyfert 2 area. One of them, J033236.72-274406.4 have IR 
and X-ray detection. 

For 4 objects, H(j, and [Nn] measurements are available and 
we have checked the log([Nii]/l 6584/H„) vs log([Om] AA4959, 
5OO7/H/3) diagnostic dia gram, see fig. [6l all of th em fall in the 
H II region delimited bv lKewlev & Dopital (l2002h . 
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Fig. 6. Diagnostic diagram for the Sample A. The dots line 
sho ws the limit of t he re gion occupied by star-forming galaxies 
from lKewlev et alj (12001.) . The dashed line represents the empir- 
ical demarcation separating star-fromation galaxies from AGN 
from lKauff^mann et all (|2003|) . 
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3.4. Estimating the metailicity of ttie ionized gas 

We have evaluated the metal abundance in the warm ionized gas. 
Such estimate reflects the current metailicity of the gas from 
which the next generation of stars will form. Accurate metai- 
licity measurements in the gas phase require the determination 
of the electron temperature Te, usually given by the ratio of auro- 
ral to nebular line flux, such as [Oiii]^^4959,5007/[Oiii]^4363. 
Unfortunately, in metal-rich environments the [Oiii]/l4363 line 
is too weak to be detected. Furthermore, even at low metallici- 
ties this auroral line is hardly detectable in low S/N spectra from 
high-z galaxies. 

To overcome this difficulty, several estimators based on 
strong emission lines are available in literature. We have used the 
R2i parameter defined by Pagel et al. ( 1979) as /?23=([0 ii].i3727 
+ [Oiii]/l/l4959,5007)/H/3. there are a large number of theo- 
retical and empirical calibrations linking R23, to o xygen abun- 
dance (jPilyugin 2001; Kewlev & Dopita 2002; Tre monti et alj 
12004 '; Pilyugin & Thuan 2005; Liang et al. 2007). Abundances 
determined by different indicators give substantial biases and 
discrepancies. For example, the difference between calibra- 
tions based on electronic temperature and on photo-ionization 
model can reach 0.5 dex jLiang et alj|2007l; iRupkeet al.1l2008t 
|KewleiL& Ellison 2008). As the aim of this work is to compare 
metallicit ies at different redshif ts, we have adopted the calibra- 
tion from lTremonti et al] (|2004|) to avoid bias related to calibra- 
tion: 



12 + logiOIH) = 9.185 - 0.313X - Q.26Ax^ - 0.321x3 



(3) 



where x - log/?23- Indeed, we want to compare sample A+B 
with the SDSS sample from which the Tremonti calibration is 
based. 

The /?23 calibration presents some limitation in abundance 
determination, namely it is double valued with 12-1- log{0/H). 
At low metailicity R23 scales with metal abundance. But from 
12 + logiOIH) - 8.3, gas cooling occurs thr ough metalHc lines 
and 7?23 decreases. The calibration from Ptremonti et alj (|2004|) 
is valid only for the upper branch of the 7^23 vs 12 -1- logO/H 
relation. We have assumed that our sample of intermediate-mass 
galaxies stands in the upper branch. In fact, galaxies in the ex- 
treme end of the lower branch are extremely rare and are asso- 
ciated with dwarf galaxies. For moderate metallicities, near the 
turning point of the relation, the uncertainties in selecting the ap- 
propriate branch is smaller than the uncertainties from sky and 
extinction. We test this hypothesis with the [Nii]/[Oii] indica- 
tor There are only 4 galaxies in our sample with [Nii] mea- 
surements: aU have log/([ Nli])//([Oil]) > -1 an d therefore 
belong to the upper-branch jKewlev & Dopitall2002f) . We think 
that, with a large level of confidence, all the objects stand in the 
upper branch of the 7^ 23 vs 12 + log{0/H). Nine objects have 
log R23 > 1 where the iTremonti et al.l (|2004|) calibration is not 
defined. In such case we have adopted the limit given where 
log7;23 = 1- 

For the 5 lowest redshift galaxies (type E3 in Table 
2), the [On] emission line falls out the wavelength range. 
We have used in this ca se the Rj, parameter d efined by 
/?3=1.35x([Oiii]^5007/Hy8) (lEdmunds & PagelllT984l) . The oxy- 
gen abundance have been estimated with the calibration pro- 
posed by Vacca & Conti ( 1992): 



12 + logiOIH) = -0.60 X log Ri - 3.24 , 



(4) 




10 11 12 

Log(M 3,e„„) 

Fig. 7. Stellar mass-metallic ity relation for the SPSS galaxies in 
contours and dots (z=0, Tre monti et alj j2004 . lBeU et al.l (l2003h 
see text ) and the 88 intermediate-z galaxies: sample A in red 
circle and sample B in red cross. Galaxies with masses under 
the completeness limit logMsteiiai < 10 are plotted with open 
symbols. The objects in sample A with log7?23 > 1 are plotted 
as small pink triangle and those detected as possible AGN are 
plotted as blue triangles. The typical error bars for bin 1, bin 2 
and bin 3 are plotted in the right. The red line indicates the me- 
dian of the intermediate-z relation.The green line represents the 
local relation fit found by 'Trem onti et al.l (|2004|) . The shift be- 
tween the local relation and our data sample is A[12H-log(0/H)] = 
-0.31dex+0.03. 



3.5. Error budget 

The uncertainties in the data were assumed to be Poisson dis- 
tributed. The error budget of emission line flux includes contri- 
butions from the sky and the object. We have propagated this er- 
ror on the Hy/Hy6 line ratio to estimate the error on Av(Balmer). 
For Av(IR) the error has been estimated taking into account the 
error on the IR luminosity and the Poisson noise of the Hj6 line. 
As the error on Av(Balmer) is dominated by the error on the Hy 
line, we can assume that the errors on Av(Balmer) and Av(IR) 
are not correlated. Thus, the error of the final extinction depends 
on the error of the two extinction estimates. The errors of line 
ratios and metailicity have been calculated by Monte-Carlo sim- 
ulation taking into account the error on extinction and the sky 
error on lines flux. 



Stellar mass 
evolution 



metailicity relation and its 



We have not found any evidence of systematic bias between the 
R2i and the Rt, calibrations. Metal abundance and error associ- 
ated are shown in Table [8] 



4.1. Tine Steilar mass - metaiiicity reiation atz ~ 0.7 

We compared the metal abundance of our sample of 88 
distant galaxies with those of local starbursts from SDSS 
(Tremonti et al. 2004). Stellar masses have been estimated us- 
ing absolute K band magnitu de with Msteiiar/^ K expected by ob- 
served rest frame B-V color dBell et al.ll2003h and converted to 
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a lSalpeteH(ll955h IM F. The convert ion between iKroupal (1200 Ih 
and ISalpeteJ d 1955 1) is described in lBell et al] (l2003l). We have 
found that distant galaxies are metal deficient compared to local 
starbursts, as shown in Fig. Q Given the small range in stellar 
mass covered by our sample, it was not possible to constrain the 
evolution of the shape of the M-Z relation at z~0.7. Thereafter, 
we have assumed that the slope of the M-Z relation has remained 
unchanged com pared to the loc al relation. At z~0.7 the local re- 
lation found bv lTremonti et al.l (j2004) is shifted to lower metal- 
hcity by A[12-Hlog(O/H)]=0.31 dex+0.03. For each galaxy we 
have calculated the diffe rence between i ts me tallicity and the 
metallicity given by the iTremonti et al.l (l2004t) relation for its 
stellar mass. The offset of the z~0.7 relation is the median of 
this difference. Only galaxies with a stellar mass over the limit 
of completeness have been used. The standard error has been 
estimated using a bootstrap method. 

The shift is assum ed to be only due to a metallicity evolution. 
Indeed, iLiang et aTl ( l2006i) argue that the evolution in the M-Z 
relation is due to a diminution of the metal content in galaxies 
rather than an increase of their stellar mass. In fact, the amount 
of star formation necessary to increase by ten times the stellar 
mass would imply that almost all the light in distant galaxies 
is associated with very young stellar population. However, the 
spectra of these galaxies show a mixture of young, intermediate 
and old stellar populations. 

There is a large dispersion of the data around the median re- 
lation at z~0.7, about +0.45 dex. To verify if this dispersion is 
intrinsic to the objects or due to the error on our data, we have 
divided our sample into three bins. Bin 2 (3) includes the 25% 
of object which are the most over (under) metal abundant com- 
pared to the median. The rest of the objects constitutes bin 1 . We 
have calculated the median error of the metallicity in each bin 
and found : 0.05, 0.03 and 0.1 (bin 1, bin 2 and bin 3). We have 
plotted the mean error of each bin in Fig. 7. There is an evolu- 
tion of the error of the metallicity through the 3 bins, which is 
well explained by the shape of the 12-i-log(0/H) vs R^i calibra- 
tion. Indeed, the metallicity varies rapidly with /?23 in metallic- 
ity range near the turn over of the function, at low metallicities. 
Nevertheless, the dispersion of the data in the M-Z relation is 
much larger that the error on metallicity. The dispersion is then 
intrinsic to the objects. We discuss about the origin of this dis- 
persion in Sect. 6.4. 

We have searched for systematic effects due to morphologi- 
cal and environment properties. First we have isolated the sub- 
sample of galaxies overstanding within in the CDFS structures. 
Metal abundance may be higher in these galaxies because of 
higher merger rate and then rapid galaxy evolution. Nevertheless 
we have not found any evidence of bias. This is probably due 
to an effect of data selection. The working sample majority in- 
cludes starburst galaxies with intermediate mass. The evolved 
galaxies in structures are expected to be massive galaxy with 
poor star-formation. 

We have also tested the effects of morphology in the locus 
of galaxies in th e M-Z relation. We have selected spiral galaxies 
in sample Ah-B (IZheng et al.ll2005b and have plotted them in the 
M-Z plot, but no effect has been found. 

Finall y, we have focused particularly on the study of LIRGs. 
Recently. iRupke et al.l (12008 !) found that LIRGs and ULIRGs are 
metal under-abundant by a factor two compared to emission-line 
galaxies in the local Universe. We have looked for an evidence 
in the sub-sample of 20 LIRGs representing 20% of the sam- 
ple A-hB. We have compared the LIRG metallicity with the star- 
burst galaxies at z~0.7 and we did not find any sh ift. The small 
metal overabundance found by iLiang et akl (|2006|) was due to a 



Table 3. Offset between local relation and the 3 redshift bin. 



Redshift bin 


A[12 + log(OIH)] 


errorA[l2 + log(0/H)] 


0.40 <z< 0.55 


-0.26 


0.08 


0.55 <z< 0.75 


-0.30 


0.04 


0.75 < z < 0.95 


-0.39 


0.09 



selection effect: LIRGs in sample B are more massive and have 
higher redshift than in sample A. Nevertheless, we have noticed 
that the z~0.7 LIRGs fall in the same locus as the local LIRGs 
in the M-Z relation. We will discuss these results in the section 
6. 



4.2. Evolution to high redshift of the stellar mass - metallicity 
relation 

One of the aims of our work is to robustly establish the evolution 
of the M-Z relation. We have first searched for a redshift evolu- 
tion of the metal content in galaxies inside the sample of inter- 
mediate redshift galaxies. As sample A+B covers a large range 
in redshift, from z=0.4 to z=0.95, we have split it into 3 redshift 
bins: 0.4 < z < 0.55, 0.55 < z < 0.75 and 0.75 < z < 0.95. 
In each redshift bin, we have measured the median shift of the 
local relation toward lower metallicity. The error of the median 
offset has been calculated by a bootstrap method. Fig. [8] shows 
the evolution of the M-Z relation along the 3 redshift bins. The 
median offsets from the local relation in the 3 bins are given in 
table[3] We find an evolution of the M-Z relation from z~0.45 to 
z~0.85 of A[12+log(O/H)]=-0.13dex. 

We have then searched for an evolution of the oxygen abun- 
dance at higher redshifts. It is important to notice that the com- 
parison between different samples needs consistent measure- 
ments of oxygen abundance and stellar mass. In fact, particular 
care has to be taken when comparing metallicities from different 
metal estimators. Especially when the aim is to evaluate the off- 
set between samples due to a evolution in metallicity. Remember 
that the systematic error due to the use of different calibration 
can range up to 0.5 dex and thu s can be higher than the metal 
evolution, cf. lRupke et al.l ( 120081) . We describe below the high-z 
samples from the hterature which have been used for this com- 
parison. 



- iLiu et al.l (l2008l) have recently published the mass- 
metallicity relation for 20 galaxies from the Deep 2 Galaxy 
Redshift Survey at 1.0 < z < 1.5. We have only considered 
the 7 galaxy spectra having individual oxygen mesurements. 
As only [O iii], li/3, [N ii] and Ha emission line measurements 
were available from the NIRSPEC/Keck observations, the 
N2 parameter has been used with the Pettini & Pag el (2004|) 
calibrator to estimate metallicity. The stellar masses were 
calculated from K-band photometry and the procedure pro- 
posed by Bundy et al. (2005X 

- At higher redshift, z~2, IShaplev et all (|2004|) observed 7 
star-forming galaxies. The oxygen abundan ce has been also 
estimated using the N2 parameter and the iPettini & Pagell 
( 12004 ) calibrator. The stellar masses have been calculated 
by fitting the spectral energy distribution to multi-band pho- 
tometry. 

- At the same redshift, (z)=2.26, lErb et all (|2006|) have mea- 
sured the metal abundance of 6 composite spectra from 87 
star-forming galaxies selected by ultraviolet rest-frame and 
binned by stellar mass. They have used the same calibration 
for 12 H- log (0/H) estimation and used the same stellar mass 
estimates as lShaplev et alj (|2004 . 
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Fig. 8. Evolution of tlie M-Z relation inside the Sample A+B. Left Panel:The 17 galaxies from the 0.4 < z < 0.55 bin. Middle 
Panel:The 38 galaxies from the 0.55 < z < 0.75 bin. Right Panel: The 15 galaxies from the 0.75 < z < 0.95 bin. The galaxies 
with masses under the completeness limit are plotted with open circles. Objects with logiRlJ) > 1 and AGN are not plotted. The 
median of each bin is marked in dashed red line, see table[3]for the A[12+log(0/H)]. The local relation from Tremonti et al. (200^ 
is marked with a green line. Only galaxies with over the completeness limit log Msteiiar < 10 have been taken to calculate the median 
shift. The two black lines in the right panel are the median of the two previous bin. 



- iMaiolino et al.l (l2008l) have recently presented preliminary 
results from the AMAZE large program which aim to con- 
strain the M-Z relation at z=3. The oxygen abundance has 
been measured by th e j^23 parameter and following the pro- 
cedure described bv 'N agao et alj (|2006|) . We have only se- 
lected the 5 objects above our stellar mass completeness limit 
and have used the 7 ^23 measurement to e stimate the metal 
abundance using the lTremonti et alj (|2004|) calibration. 

Three samples use the N2 parameter to evaluate the oxy- 
gen abundance. Unfortunately, metallicity estimation based on 
N2 and 7?23 can give strongly discrepant results and thus can- 
not be directly compared. The N2 parameter has several disad- 
vantages: the [Nii]/H„ ratio saturates at 12 4- log(0/H) > 8.8 
and becomes insensitive the metal abundance. Furthermore, it is 
very sensitive to the ionization parameter and the contribution 
of the primary and second nitrogen whose evolution is not yet 
well constrained. To compare the galaxy metallicities using the 
N2 para meter we hav e decided to re-evaluate the metallicities 
using the lDenicolo et a l. (2002) calib ration . Indeed, this calibra- 
tion is similar to the llYemonti et^ (|2004|) and other R23 based 
calibration in the working range 8.4 < 12 4- log(0/H) > 8.7 
dKewlev & EUisonll2008tlTremonti et"ani20 04V For the galaxies 
in the sample A-l-B having [Nn] measurements, we have com- 
pared the metallici ty to those of several N 2 calibrations. We have 
confirmed that the lDenicolo et alj (|2002|) calibration give similar 
results than the abundances estimated in this work. The M-Z di- 
agram for the local, intermediate and high redshift samples are 
plotted in Fig. 9. 

We have measured the median offset of the local relation and 
the associated error in all the high-z samples. The offset A(0/H) 
for the 3 redshift bins of the intermediate galaxies and for the 
four high-z samples are plotted in Fig. [TO]as a function of the 
lookback time. We found that the evolution of the mean metal- 



licities from local Universe to a lookback time of ~12Gyrs is 
linear and with a slope given by: 



logZtb/Zo = -0.046 xfG 



yrs : 



(5) 



If the iPettini & Pagell ( |2004|) is used for metalUcity estimate of 
high-z samples, the evolution remains linear with a slope of 
-0.051. For comparis on we o verplot the metallicity evolution 
extrapolated bv ISavaglio et al.l (12005). The empirical evolution 
model has been constructed by shifting the local M-Z relation 
in log Msteiiar Contrary to current work where the shift has been 
evaluated in term of metallicity for the reason explained above. 
Their empirical model is concordant with a close box model with 
an exponential star formation efficiency. We discuss in sections 
5 and 6 the di s crepan cy between our observations and those of 
ISavagUo et alj (I2005h . as well as the consequences for our un- 
derstanding of the evolution of the M-Z relation. 



5. Possible systematic uncertainties 

In this section we discuss the possible systematic errors that can 
affect our determination of the evolution of the M-Z relation. 



5.1. Metallicity calibrators 

As we have emphasized in §3, the determination of metallic- 
ity using dif ferent calibrations ma y give very discrepant results. 
However, IKewlev & EllisonI (i2008l) have shown that the relative 
metallicities between galaxies can be reliable estimated within 
~0.15dex when using the same strong line calibration. In this 
work we have taken a special care to compare our sample to the 
local g alaxies us ing the same metallicity calibration and method- 
ology. |LiuetaL| (|^0l) have compared the M-Z relation at z~l 
to the local relation using N2 for both redshifts. They have found 
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Fig. 9. Stellar mass-met allicity relatio n for the SDSS galaxies in 
contours and dots (z=0, jTremonti et a l. 2004; Bell et al. 2003) ) 
and the 58 intermediate-z galaxies above the completeness mass 
limit in red fill circle. Galaxies under the completeness mass 
limit a re plotted in red open circle. The local relation is plotted in 
green jTremonti et al.ll2004b and the median of the inte rmediate- 
z relat ion is in red. The 7 individual galaxies at z>l of iLiu et all 
(1200 8') are plotted as black triangle. The two z>2 samples are 
plotted in b lue: inverted b l ue tria ngles indicate the 7 star-forming 
galaxies of lShaplev et al. ( 20041) : blue stars indicate the 5 stellar 
mass bin of lErb et alj (l2006l) . The z > 3 sample of Maiolino et 
al. (2008) are plotted in green squares. 



that z~l galaxies have metal abundance 0.22 dex lower than in 
present Universe. This is quite comparable to our results. 

The comparison between the metallicity at z~0.65 and 
higher redshift samples is less robust: low S/N spectra and differ- 
ent metal indicators have been used. In fact, high-z samples use 
usually the N2 indicator. Several authors have tried to build con- 
versi ons between metallicity calibrations, e.g. lKewlev & EUisonl 
(12008.) . However, these conversions have been derive d from lo- 
cal sam ples and may be not valid at higher redshift. iLiu et alj 
(|2008|) have suggested that at z~l the metallicities of galaxies 
are overestimated by ~0.16dex due to different physical con- 
ditions in HII regions at higher redshift: violent star-formation, 
shocks, higher ionization parameter. No evidence of such dif- 
ferent conditions of HII regions have been found in our sample, 
according to the diagnostic diagrams. 



5.2. Aperture effect 

The eff'ect of the aperture when comparing intermediate red- 
shift data to the SPSS sample has b e en widely studied 
by several authors (ISavaelio et all 12 0051: Liang et al.l l2004at 
iKobulnickv & KewlevI 120041: iKewle v & Ellison 2008). The 
SDSS observed only the center part of galaxies with a 3" 
aperture fiber. The metallicities estimated are expected to 
be higher than the average metallicity of the galaxies be- 
cause of the metallicity gradient with the distance from 



-0.5 



-0.6 



"> 1 r- 



Savaglio et al. 



[] 







10 



12 



2 4 6 8 

Lookback time (Gyrs) 

Fig. 10. The metallicity shift from the local relation of the 4 
high-z sample as a function of the lookback time is plotted as 
black open squares. The metallicity shift for the 3 redshift bin 
of the IMAGES sample are in black squares and the median of 
the 3 bins as a red triangle. The red line is the linear fit of the 
5 high-z data poin ts verifying Ari2-Hlog (O/H)]=0 at z=0. The 
evolution found bv lSavagho et alj (l2005l) in the frame of a close 
box model is indicated with a dashed line for log Msteiiar= 10.45 
galaxies. 



the cen ter. dyila-Costas & Edmundsl ll992l: IZaritskv et all 1 1994 
TGarnett etalJ 119971: iRoUeston et all l2000h . iKewlev & Elhsonl 



(i2008) have shown that the effect is relevant for massive galax- 
ies Msteiiar > IO'^Mq. Howcvcr, the bias due to aperture ef fect in 
the SDSS is not larger than Z~ 0.1 dex dKewlev et al.ll200 5) and 
thus have a minor impact on the estimation of the evolution of 
the M-Z relation. 



5.3. Extinction 

The effect of extinction on the estimation of metallicity is one 
of the main systematic errors affecting the M-Z relation. For 
example, an underestimation of 1.5 mag in extinction leads to 
an overestimation of ~0.15dex in oxygen abundance. Several 
works have assumed an average extinction to correct spectra of 
galaxy delimiting a large range in stellar m ass. However, ther e 
is a correlation between mass and extinction dLiang et al.ll2007h : 
massive galaxies have higher extinction. Assuming a median ex- 
tinction leads to an overestimation of the metallicity in massive 
galaxies and an underestimation in low mass system. At z~0.6, 
LIRGs which have high content of dust and gas represent 15% 
of intermediate mass galaxy population. Extinction is expected 
to have a large impact on any abundance estimate at this redshift. 

The effect of extinction may be one of the reasons for the 
differe nce between our observations and those of ISavaglio et al.l 
(2005"). They have observed that the average metallicity of mas- 
sive galaxies has not evolved from z~0.6 to present days and thus 
have claimed evidence for downsizing in metalhcity. However , 
the majority of massive galaxies in the ISavagho et alj (l2005h 
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sample come from measurements of iLillv et alj (l2003h which 
are from spectra with low S/N and spectral resolution. Moreover, 
they have applied an average Ay = 1 to correct emission lines. 
This extinction value is underestimated b y about ~lmag for 
at least one third o f the Lilly et al.* (2003) sample, which cor- 
responds to LIRGs. [Liang et arr (2006) h ave estimated th at the 
underestimation of the extinction in the L illy et al] ( 12003 ) sam- 
ple leads to an overestimation of ~0.3 dex on metallicity, which 
can explain the discr epancy with our sam ple. The galaxies in 
the GDDS sample in ISavaglio et alj (l2005i) have low luminosi- 
ties and can be also affected by the low S/N of the spectra 
which does not allow to retrieve proper extinction from Balmer 
lines for individual galaxies. These effects are amplified for 
evolved massive galaxies especially those experiencing succes- 
sive bursts and containing a substantial fraction of A and F stars. 
Unfortunately these stars are predominant in intermediate mass 
galaxies, and thus contaminate extinction measurement by their 
large B almer absorption. Fur thermore, the average extinction 
used by ISavaglio et al.l (l2005h can be underestimated because 
the combining spectra used to measure Hy and H/3 lines is dom- 
inated by low extinction systems. 

5.4. Stellar mass 

A bias in the estimation of stellar mass at high redshift must 
be considered. In this work, stellar mass es have been esti mated 
from Msteiiai / I^K ratio using the method of lBell et"an ( l2003l) . This 
method uses the tigh correlation between rest frame optical color 
and Msteiiai /i'K ratios. The amount of light due to red giant stars 
is corrected in the Msteiiar/^K using the g-r colors. At higher red- 
shift, the influence of TP-AGB stars in the derivation of the stel- 
lar masses could result in an overestimation of th e stellar mass by 
-0.14 dex jMaraston et alj|2006l: [Pozzetti et al.ll2007l) . However 
this bias has only a thin effect on the evolution of the M-Z rela- 
tion. The bias in stellar mass directly translated into a systematic 
effect of the evolution of the M-Z relation is around -0.05 dex. 



5.5. Uncertainties budget 

We summarize all possible sources of systematic uncertainty 
identified so far in Tables 4 and 5. We expressed all of them in 
terms of their influence in the shift of the M-Z relation from local 
Universe to higher redshift. In the intermediate redshift sample 
we have taken a particular care to diminish the bias from extinc- 
tion estimates and from metallicity calibrations. In the first place, 
the methodology has allowed us to reduce the bias due to ex- 
tinction. Secondly, we have used a similar methodology and the 
sa me metallicity c alibrators for our sample and the local sample 
of iTremonti et al.l (2004). The remaining systematics may come 
from stellar mass estimates and from the aperture effect on the 
SDSS data. However, as illustrated in Table 4, the systematic er- 
ror from these two sources is small and thus does not change our 
results. 

The shift of the M-Z relation from intermediate to high red- 
shift is less secure. The sources of systematic error come from 
the use of different metallicity calibrations and a possible evolu- 
tion of the HII properties at high-z, see Table 5. We have tried 
to minimize the bias from metallicity estimates selecting for the 
high-z sample, the metallicity calibration that b etter follows the 
one us ed in this work. However, even using the Pettini & Pagel' 
(|2004|) calibration instead of the Denicolo et al. (2002) calibra- 
tion, we do find a linear evolution of the metal content in galax- 
ies and with a similar slope. As the evolution of the M-Z relation 



Table 4. Identified systematic uncertainties that could impact on 
the shift of the M-Z relation between local and z~0.6 galaxies. 



Possible bias 


Z 


Comments 


Aperture effect 


0.1 


See Sect. 5.2 


Stellar mass 


-0.05 


See Sect. 5.4 



is strongly constrained at z~0.6, the bias in the high-z sample 
will not change the shape of the evolution function but only the 
slope. Even accounting for all the systematics shown in Table 
5. the observations from intermediate and high redshift samples 
are no t compatible with the evolution found by .Savaglio et al] 
(120051) . 



6. Discussion 

6.1. Discarding the closed-box model 

We first compare our observations w ith a simple model of 
galaxy evolut ion: the closed-box model (ISearle & Sargent! 19721; 
iTinslevI [l980). In this model galaxies are isolated systems 
which evolve passively transforming their gas into stars, en- 
riching their content in metals. The amount of gas converted 
into stars from z~0.6 to present can be est imated, following 
the same methodology as iLiang et al] ( l2006l) and proposed by 
iKobulnicky & Kewleyl (l2004l) . The metal abundance Z is related 



to the gas mass fraction by 
Z(f) = Y \nl/fi{t). 



(6) 



where jj = Mgas/(Mgas + Msteiiai) is the fraction of gas. As the 
yield is constant the variation of metallicity depends only on the 
variation of the gas fraction : 



d(logZ)/dyu = 0.434/yuln|/, 



(7) 



We assume that local galaxies have a gas fraction /i()=10%. The 
metallicity of the gas phase has decrease by two from z=0 to 
z=0.6. Then by equation (7) the gas fraction at z=0.6 predicted 
by the close box model is around mu=30%. This means that the 
present day intermediate galaxies would have converted 20% 
of their present b aryonic mass, from gas into stars. Recently, 
iPuech etal] (12008*) have studied the evolution of the Tully-Fisher 
relation at z~0.6 using a representative sample of 65 emissions 
line intermediate galaxies. They have found a significant evolu- 
tion of the K-band TFR zero point, which they attribute to an 
average brightening by 0.66 mag. After considering several al- 
ternatives for the evolution of the stellar mass, they conclude 
that on average, spiral galaxies have doubled their stellar masses 
since z=0.6. The closed-box model predicts a much lower effi- 
ciency in converting gas into stars during the same elapsed time. 
There would be not enough gas from the close box model to feed 
the stellar mass doubling. This discrepancy can be explained if 
~30% of the stellar mas s was formed by external gas supply 
probably metal enriched dEdmundsll 19901) . We can definitively 
rule out a closed-box model as a viable description of the metal 
content in galaxies from z=0.6 to z=0. 

The closed-box model has been taken as a first approach 
to exp lain the metal evolution of galaxies, e.g lSavaglio et al.l 
(120051) . However, we notice that the close box model is not 
consistent with the two major scenarios of galaxy evolution. 
Indeed, to resolve the G-dwarf problem and to explain the evo- 
lution of disk, the secular scenario needs to take into account 
large amount of infalling gas, e.g. primordial gas in filaments 



12 



M. Rodrigues et al.: IMAGES IV: Strong evolution of the oxygen abundance 



dSemelin & Combesll2005h . In the hierarchical scenario the evo- Thus existing modeling have difficulties in reproducing our 
lution of galaxies is driven by mergers and interactions, thus by observations, i.e. a strong evolution of the metal abundance of 
gas exchanges. intermediate mass galaxies over the last 6 Gyrs. 



6.2. Comparison with previous models of galaxy chemical 
evolution 

We compare our results with galax y evol ut ion simulations 
by Ide Rossi et al.l (l2007l) : iBuat et al.1 (|2008|) : iMouhcine et"aLl 
(12008'). We have only compared the evolution of metal content 
in the gas phase founded by the models in the same mass range 
than in our sample. 

Within the fra mework of the hierarchical scenario, 
Ide Rossi et al.l (l2007h have simulated the evolution of the M-Z 
relation from present to z=3 . They predict very little evolution on 
the relation of only -0.05 dex since z=3 and that massive galax- 
ies with Msteiiar > 3 X IO^'Mq havc reached the local relation 
well before z-l. Thus they predict that the spectra of massive 
galaxies should be mainly composed of an old metal-rich stellar 
population. This is in strong contradiction with our results. In our 
study 38% of the galaxies are above the critical mass and have 
strongly increased their metal content, whic h is consistent wit h 
the evolution of the Tully Fisher relation dPuech et al.l l2008 l). 
Moreover, the results of stellar continuum fitting indicates that 
the predominant stellar population is composed of A and F stars, 
which suggest a significant stellar formation 1 .5 Gyr prior to our 
observations. The predominance of intermediate age stellar pop- 
ulation in galaxies at intermediate redshifts had been previously 
poi nted ou t by Hammer et al. (2001). 

iMouhcine et alj (l2008l) have simulated numerical models of 
galaxy evolution in the hierarchical scenario and they also pre- 
dict an important downsizing effect. Metal abundances of galax- 
ies with Msteiiai > IO'^Mq do not evolve since z~1.2. However, 
they have found a stellar mass assembly and a metallicity disper- 
sion consistent with our results. 

Finally we compared our observations wit h the chemical 
evolution model presented bv iBuat et alj ( l2008l) . They investi- 
gate the star formation history of star-forming galaxies from z=0 
to z=l and compare with evolution models developed for local 
spiral galaxies within the frame of a secular evolutionary sce- 
nario. From a model they predict the SFR and the M-Z relation 
evolution for intermediate mass galaxies. Their predicted evolu- 
tion of metallicity is half of our observations: they predict a shift 

of 0.2 dex at z=l (against 0.4 dex in this study, see Table 

3). 

They also predict that at z~l massive galaxies have in- 
creased more rapidly their metal content than low mass galax- 
ies wh ich is at odds with most current models. The iBuat et alj 
(|2008|) model is based on the evolutionary history of the Milky 
Way, a galaxy which may not be representative of typical spirals 
dHammer et al.ll2007l) . For example, it is unlikely that the Milky 
Way has experimented in the past 8 Gyrs a strong episode of star- 
formation such as a LIRG, which is quite common in galaxies of 
our sample. 

At z= 0.4-0.9, the fraction o f emission line galaxies is 60- 
70% ("e.g. lHammer et afl dlQQTl) ). As such, any deviations from 
local values in a sample of emission line galaxies is at odd with 
the above modelling. If we constrain Sample A to the redshift 
range from 0.4 to 0.7 (limit technical problems like lines out of 
range), we find that 96/136 galaxies are emission lines galaxies 
(31 withEQW(OII)<15and65 with EQW(OII)>15) and only 40 
are absorption line galaxies. From the 65 emission lines galaxies 
with EQW(OII)> 15, 40 have been studied the rest being rejected 
mainly for instrumental reasons (sky lines, instrument fakes etc). 



6.3. The star-formation efficiency of intermediate mass 
galaxies atz~0.6 

We estimate the doubling time, defined by - Msteiim/'S'^^^iR, 
for emission line galaxies in the IMAGES sample. The SFR for 
galaxies without IR detection have been estimated using the Lir 
limit of detection of MlPS/Spitzer and the 5'F/?28ooa- To de- 
termine the unobscured S ^^^28{)oa' \i&sd the UV calibration 
ofl Kennicutddl998l) and the rest-frame 2800A luminosity from 
SED fitting. 

We confirm that massive galaxies convert less gas into stars 
than lower mass galaxies at 0.5 < z < 1: e.g the to for a mas- 
sive galaxy log Msteiiai = l 1 is above 3 Gyrs whereas a less mas- 
sive galaxy logMsteiiai = 10.3 doubles its stellar content in less 
of 1 Gyrs. This effect is usually called downsizing: low mass 
galaxies are forming stars at latter epoch than massive galax- 
ies. Severals authors have pointed out that downsizing is in con- 
tradiction with the hierarchical scenario. Massive galaxies form 
later than low mass galaxies by accretion of small dark halos, in 
the bottom-up scenario. However Mouri & Taniguchi (2006) and 
iNeistein et al.l d2006i) have demonstrated that the downsizing of 
star-forming galaxies is inherent in gravitational processes of the 
hierarchical scenario. Low mass galaxies have very long gravita- 
tional collapse time-scales due to their small potential well. Thus 
they start to produce stars later than massive galaxies. 

The small values of doubling time in intermediate mass 
galaxies at z~0.6 (see Figure 12) confirm the strong evolution 
of these sources at this epoch. We do find that the doubling time 
diminishes towards redshift which is due to the star-formation 
increase with redsh ift.These obs ervations are in good agreement 
with the results of iPuech et al.l (|2008) who argue that galaxies 
have doubled their stellar mass in the last 6 Gyrs. 



6.4. Similarities between local LIRGs and distant starbursts 
& LIRGs 

The IMAGES sample has enabled us to investigate several prop- 
erties of distant starbursts, such as metallicity and star-formation 
rates, etc. These observations reveal that local LIRGs and distant 
starbursts/LIRGs share many properties: 

- Same locu s in the M-Z relation 

iRupke et~a l. (2008) have shown that local LIRGs do not fol- 
low the local M-Z relation. Recall that the average SDSS 
galaxies populating the local M-Z plane are forming stars at 
moderate rate (less than few solar masses per year), much 
lower than for starbursts and LIRGs. Local LIRGs (see 
Figure 11) have a deficiency in metal compared to other lo- 
cal galaxies and delineate the same locus in the M-Z relation 
than the high-z starbursts & LIRGs. Both populations show 
a large and similar dispersion of their metallicity at a given 
stellar mass. Thus local LIRGs and distant starbursts share a 
same metal content at a given stellar mass. Notice that in this 
work, we have not found differences between distant LIRGs 
and distant starbursts. 

- Similar high star-formation efficiency 

Fig. 12 shows the distribution of doubling times for local 
LIRGs (Rupke et al. 2008) and IMAGES galaxies. In order 
to properly compare the two samples, we used a same stellar 
mass range to avoid biases (e.g. lower doubling time from 
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Fig. 11. M-Z rel ation for starburst galaxies at z~0.6 and local 
LIRGs from Rup ke et all (jToOS). The symbols used in this plot 
are the same as in fig. 9. The local LIRGs are plotted as black 
stars. 



low mass galaxies in the intermediate mass galaxy sample). 
The two distributions are very similar. The Kolmogorov- 
Smirnov test attests that the probability that the two distri- 
butions arise from the same population is about 40%. We 
notice that stellar mass for local LIRGs has been estimated 
from their absolute magnitude Mk and assuming a MjL^- 1 . 
This crude assumption can lead to a slight overestimation of 
their stellar mass (by ~0. 1 dex) although it does not affect the 
overall result. 

- Similar morphologies: dominated by mergers and spi- 
rals 

Severals studies in the local Universe have pointed 
out the link betw ee n LI RG s, merging of g as-rich 
spirals! VeilleuxetiD 120021: ISanders & Ishidal 120041; 
IWang et al.,,20^6l) . At higher redshifts a significant fraction 
of the dista nt galaxies have morphologies that inclu de spiral 
to peculiar ("Zheng et al]|2005l: iNeichel et alll2008l) . Distant 
LIRGs are also associated with merger at different sta ges or 
with large disks dZheng et alJl2005l;lFlores e"tanil999h . 

Fundamental properties such as morphologies, star-formation 
efficiency and mass-metallicity are similar for both local LIRGs 
and distant starburst/LIRGs. We may thus wonder whether local 
LIRGs and distant starburst also share the same physical mecha- 
nisms explaining their properties. The only discrepancy between 
the two populations is their numerical abundance: while local 
LIRGs only represent 0.5% of the local massive galaxies, dis- 
tant starbursts studied here represent approximatively 60% of 
distant massive galaxies dHammer et alj|1997l) . We suggest that 
local LIRGs are the latecomers of the main population of mas- 
sive starburst galaxies visible at z~0.6. 

What is the physical process explaining the properties of 
these two galaxy populations of actively star forming galaxies? 
Since star formation increases the abundance of a syst em, the 
only processes that can reduce it are gas motions (lEdmundsl 
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Fig. 12. His togram of the doubling time of local LIRGs from 
Rup ke et akl (|2b08) (upper panel) and emission lines galaxies 
from the IMAGES sample (botton panel). The distribution of 
IMAGES sample for galaxies delimiting the same logMsteiiar 
than local LIRGS is overplotted in red. 



1990: Koppen & Edmunds' HH iDalcantonI IIoOTI) . For local 
LIRGs, Rupk e et al.. (|2008) have proposed that the encounter be- 
tween two gas rich galaxies causes lower metal abundance gas 
from the outer regions to fall into the merger central region, pro- 
viding dilution of the metal abundance. The input of gas pow- 
ers the star formation and dilutes metals at the same time, ex- 
plai ning their l ocatio n in the M-Z relationship. Independently, 
Ha mmer et al. I (l200l have suggested a scenario in which present 
spiral galaxies are built by the gradual infall of the gas being ex- 
pelled during the first stages of a major merger. Both scenarios 
are based on supply of unprocessed gas from the outskirts of en- 
counters during a merger event. For both cases, outflows are un- 
likely the cause of the gas dilution as they may be dominated by 
a neutral phase with the metal abundance of the central regions 
(lRupkeetal.ll2008ll2005 t). 

Because distant starbursts & LIRGs are a common phe- 
nomenon in the distant Universe, a merger scenario to explain 
their properties may have a considerable impact for the forma- 
tion o f present-day spirals. It has prompted by iHammer et al.l 
('2005') to propose a scenario of disk rebui lding for the majo rity 
of spirals galaxies. Recent simulations of iLotz et al.l ( l2008l) in- 
dicate that the majority of mergers can reconstruct a low mass 
star-forming disk. Moreover, the large dispersion in metallicity 
we observed at intermediate redshifts is qualitatively consistent 
with the diversity of galaxy building histories within the frame- 
work of a merger scenario. 

The external gas supply through merger process is in good 
agreement with Damped Lyman Absorber (DLA) observations. 
DLAs have a lower metal content compared to emission selected 
galaxies and have a wide dispersion in metallicity. Absorption 
measurements through DLAs trace outer regions of galaxies and 
probably e xtend envelopes of unprocessed gas on the outskirts 
of galaxies. IWolfe et al.l(l2003h have proposed a scenario for the 
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formation of spiral galaxies from DLAs in which disk galax- 
ies are surrounded by an halo of neutral gas during t heir as- 
sembly . The simulation of Navarro (2004) and Govern ato et alj 
(I2OO7I) have reproduced many properties of local spiral galaxies 
by merger of gas rich encounters. 



7. Conclusions 

We have observed 88 intermediate mass galaxies with high qual- 
ity spectra and moderate resolution to establish a robust study of 
the stellar mass metallicity relation at z~0.6. Metal abundance 
have been estimated following a careful analysis of each individ- 
ual spectra: evaluation of underlying absorption in Balmer lines 
by stellar continuum subtraction, robust extinction measurement 
and evaluation of AGN contamination. 

We confirmed the shift about ~0.3 dex to lower a bundance of 
the M -Z relation at z~0.6 which has been found by ^ Liang et al] 
(120061) . Combined with the robust selection of our sample, the 
relatively high number of studied galaxies in this study allows to 
derive the evolution of [12-i-log(0/H)] over the redshiftbin 0.4 < 
z < 0.9. It evidences that the median metal abundance of galaxies 
diminishes with the redshift. Using similar measurements from 
samples at higher redshift we recover the evolution in metals 
in the ionised gas from z=0 to z~3: we find a linear evolution 
of metallicity as a function of the lookback time. This simply 
means that the evolution of the gas phase in massive galaxies is 
still active down to z=0.4 conversely to the popular belief that 
all massive galaxies have their stellar content locked since z=l. 

Our result s are th us in strong disagreement with those of 
ISavaglio et al.l (|2005|) et al. They indeed found that the evolu- 
tion of galaxies can be at first order explained by a close box 
model and that galaxies with logMsteiiar > 10'"Mo have already 
reached local metallicities at z~0.6. On contrary we have dis- 
carded the close box model as a valid approximation of galaxy 
evolution. We suggest t hat this disc repancy is due to the low S/N 
quahty of the| Savaglio et al.l (l2005h data. We notice however that 
ISavaglio et al.l (i2005) results have been widely used to construct 
metal content evolution simulations, explaining why our results 
are not reproduced by these models. We are calling for a new 
generation of modelling able to account for our observations. 

We have put into evidence the need for a significant con- 
tribution of external gas supply in the metal content evolution 
of galaxies. Since 6 Gyrs ago, a significant fraction of massive 
galaxies have undergone a rapid growth in which 30% of their 
present stellar mass has been formed from external gas supply. 
The addition of gas in galaxies explain the metal deficiency of 
star-forming galaxies at z~0.6 at all mass scales. The origin of 
the gas supply is a key question. In the secular scenario, the in- 
falling gas come from pristine gas in the large scale filament 
structures. The perturbed kinematics and morphology of the in- 
termediate redshift galaxies have lead us to favor the hierarchical 
scenario where the unprocessed gas could arise from large en- 
velopes of gas lying in the halo of galaxies after merging events. 

We have pointed out that local LIRGs and distant star- 
bursts and LIRGs share similar properties such as star-formation 
efficiency, metal abundance and morphology. We suggest an 
unique evolution scenario for both populations, which is linked 
to the gas infall provided by galaxy merging and interactions. 
Progenitors of spiral galaxies have their disks rebuilt through an 
infall of the gas which has been previously expelled in the first 
stage of a gas rich merger Local LIRGs are in this case latecom- 
ers of the starbursts/LIRGs population predominant at z~0.6. A 
better understanding of the kinematics and metallicity of galac- 



tic halos in gas-rich system at intermediate redshift is needed to 
disentangle the origin of the external gas supply. 
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Table 5. Identified systematic uncertainties that could impact on the shift of the M-Z relation between z~0.6 and high-z galaxies. 



Possible bias 


Z 


Comments 


Metallicity calibration 


0.3 


Systhematic between N2 and R^j calibrations. 


Properties of H n region 


-0.16 


For high-z sample with N2 calibrations. 


Extinction 


-0.3 


i?23 estimated sample without extinction correction 


Stellar mass 


-0.05 


See Sect 5.4 
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Table 6. Basic data from galaxies in sample A. The IR luminosity of galaxies has been estimated from the mid-IR catalogue of iLe Floc'h et alJ 
( 12001) and using the procedure of ICharv & ElbazI ( 1200 ih . The rest frame magnitude in J,K-band and the luminosity at 2800A have been derived 
by modeling galaxy SEDs using ISAAC and ACS multi-band photometry. 



Name 



RA (J2000) DEC (J2000) /ab 



J033210.92- 
J033211.70- 
J033212.30- 
J033212.39- 
J033212.51- 
J033213.76- 
J033214.48- 
J033215.36- 
J033217.36- 
J033217.75- 
J033219.32- 
J033219.96- 
J033222.13- 
J033223.06- 
J033223.40- 
J033224.60- 
J033225.26- 
J033225.46- 
J033225.77- 
J033226.21- 
J033227.36- 
J033227.93- 
J033227.93- 
J033229.32- 
J033229.64- 
J033229.71- 
J033230.07- 
J033230.57- 
J033231.58- 
J033232.13- 
J033232.32- 
J033232.58- 
J033233.00- 
J033233.82- 
J033233.90- 
J033234.04- 
J033234.88- 
J033234.91- 
J033236.37- 
J033236.52- 
J033236.72- 
J033236.74- 
J033237.26- 
J033237.49- 
J033237.96- 
J033238.77- 
J033238.97- 
J033240.04- 
J033240.32- 
J033243.96- 
J033244.44- 
J033245.ll- 
J033245.51- 
J033245.63- 
J033245.78- 
J033248.84- 
J033249.58- 
J033252.85- 



274722.8 
274507.6 
274513.1 
274353.6 
274454.8 
274616.6 
274320.1 
274506.9 
274307.3 
274547.7 
274514.0 
274449.8 
274344.5 
274226.3 
274316.6 
274428.1 
274524.0 
275154.6 
274459.3 
274426.3 
275015.9 
274353.6 
275235.6 
275155.4 
274242.6 
274507.2 
274534.2 
274518.2 
274612.7 
275105.5 
274343.6 
275053.9 
275030.2 
274410.0 
274237.9 
275009.7 
274440.6 
274501.9 
274543.3 
275006.4 
274406.4 
275206.9 
274610.3 
275216.1 
274652.0 
274732.1 
274630.2 
274418.6 
274722.8 
274503.5 
274819.0 
274724.0 
275031.0 
275133.0 
274812.9 
274531.5 
275203.1 
275207.9 



0.417 
0.677 
0.645 
0.422 
0.732 
0.679 
0.546 
0.860 
0.647 
0.734 
0.725 
0.784 
0.541 
0.734 
0.616 
0.538 
0.666 
0.672 
0.833 
0.495 
0.769 
0.458 
0.383 
0.510 
0.667 
0.737 
0.648 
0.679 
0.654 
0.682 
0.534 
0.670 
0.669 
0.666 
0.619 
0.703 
0.677 
0.665 
0.435 
0.689 
0.666 
0.784 
0.736 
0.423 
0.620 
0.458 
0.420 
0.523 
0.619 
0.533 
0.416 
0.436 
0.562 
0.858 
0.534 
0.278 
0.415 
0.684 



53.04549408 

53.04874039 

53.05125809 

53.05161285 

53.05210876 

53.05731201 

53.06032181 

53.06398392 

53.07234955 

53.07396698 

53.0804863 

53.08317947 

53.09218979 

53.09606552 

53.09751511 

53.10250473 

53.10525131 

53.10609818 

53.10738373 

53.10919571 

53.11399078 

53.11636734 

53.11637878 

53.12216568 

53.12350845 

53.12377548 

53.12527084 

53.12736893 

53.13157272 

53.13389587 

53.1346817 

53.13574982 

53.13750839 

53.1409111 

53.14123917 

53.14182281 

53.14534378 

53.14544678 

53.15153885 

53.15215683 

53.1529808 

53.15309906 

53.15524292 

53.15618896 

53.15816498 

53.16156006 

53.16236496 

53.16683197 

53.16801453 

53.1831665 

53.18515015 

53.18795013 

53.18963242 

53.19010544 

53.19076157 

53.20350266 

53.20658112 

53.22020721 



-27.78966522 
-27.75211334 
-27.75362968 
-27.73155403 
-27.74856567 
-27.77127075 
-27.72223854 
-27.75192451 
-27.71868706 
-27.76324844 
-27.75389862 
-27.74717331 
-27.72902679 
-27.70730209 
-27.7212677 
-27.74114418 
-27.75665855 
-27.86517906 
-27.74981499 
-27.74064445 
-27.83775711 
-27.73156548 
-27.87656403 
-27.86540031 
-27.71182251 
-27.75200462 
-27.75950241 
-27.75506783 
-27.77020454 
-27.85154152 
-27.72879028 
-27.84830666 
-27.84171295 
-27.73611832 
-27.71053696 
-27.83602524 
-27.74459839 
-27.75053024 
-27.76201439 
-27.83511353 
-27.73512459 
-27.86857033 
-27.76953125 
-27.87112999 
-27.78109741 
-27.79225922 
-27.77506256 
-27.73850822 
-27.78967285 
-27.75096512 
-27.80527496 
-27.78999901 
-27.84193993 
-27.85918045 
-27.80357933 
-27.75874901 
-27.86752319 
-27.86885071 



20.27 
22.63 
21.60 
21.40 
22.81 
22.92 
22.84 
22.29 
21.30 
22.14 
22.26 
22.60 
23.03 
22.17 
20.84 
22.05 
21.40 
21.09 
22.86 
23.28 
21.83 
23.53 
20.55 
22.25 
20.90 
22.39 
21.71 
20.81 
22.44 
22.34 
22.84 
21.86 
21.43 
21.51 
21.22 
22.35 
23.13 
22.48 
22.14 
21.65 
22.07 
22.30 
22.38 
21.09 
21.97 
21.13 
21.07 
20.89 
23.00 
22.41 
20.57 
20.81 
22.70 
22.40 
21.70 
21.50 
21.14 
23.06 



10.95 

10.90 
10.94 



11.35 - 

11.27 
11.16 

11.01 

10.69 

11.77 

11.41 
11.07 

11.62 

10.92 

11.21 

10.88 
10.98 
11.05 
11.12 

10.47 - 

11.48 
10.35 
10.89 

10.79 - 



9.81 

10.70 

11.10 



21.97 
19.89 
21.89 
21.59 
19.87 
19.41 
19.28 
21.57 
21.77 
21.20 
21.25 
20.88 
19.03 
21.43 
23.00 
20.47 
21.63 
22.68 
20.94 
18.59 
21.98 
17.86 
21.59 
19.42 
22.74 
20.93 
21.44 
23.00 
20.35 
21.71 
19.26 
21.54 
23.18 
22.61 
21.84 
20.69 
20.65 
20.42 
21.04 
21.78 
22.02 
21.39 
21.15 
20.87 
20.53 
21.32 
21.22 
22.02 
19.62 
20.16 
22.02 
22.08 
19.52 
21.02 
21.56 
19.56 
21.48 
20.64 



-21.89 
-19.74 
-21.84 
-21.57 
-19.56 
-19.07 
-19.17 
-21.47 
-21.67 
-21.11 
-21.16 
-20.78 
-18.82 
-21.33 
-22.98 
-20.38 
-21.54 
-22.60 
-20.71 
-18.38 
-21.89 
-17.59 
-21.50 
-19.09 
-22.66 
-20.82 
-21.34 
-22.91 
-20.26 
-21.68 
-19.15 
-21.46 
-23.15 
-22.55 
-21.74 
-20.60 
-20.56 
-20.31 
-21.01 
-21.69 
-21.94 
-21.30 
-21.06 
-20.78 
-20.43 
-21.24 
-21.13 
-21.93 
-19.53 
-20.07 
-21.95 
-22.05 
-19.43 
-20.68 
-21.51 
-19.47 
-21.45 
-20.57 
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Table 7. The extinction Av(Balmer), Av(IR) and the adopted Ay. The Av(lR) calculated with the lower limit of IR detection are upper limit. The 
final extinction estimated only with the limit Ljr are noted with '. Three SFR are given : SFRi^, SFRjit,, 5Fi?uv- The SFR^n, is corrected from 
extinction and aperture. To determine the unobscured 5F/fuv, we used the rest-frame 2800A luminosity from SED fitting 



Name 



Ay 

B aimer 



IR 



SFRjR 



S F/?Hb 

Meyr-l 



S FRuY 
Moyrl 



-Q-ff+UTT 



1 02+"-^'* 

^ -^^-0.25 

0.31' 

1 57+°-35 

^•-"-0.36 

1.97' 

1.55 + 0.26 
0.28^ 
1.70' 

2 16+°-47 



11.08-- 

< 16.58 

i8.9i!;;j:i 

-i 4Q+4.89 

< 16.65 

59.06!ii^ 

< 14.59 

< 20.79 

47 Qc+37.68 
^' ■^-'-40.91 

< 25.25 

o 1 Q+5.06 
3 13 

< 20.77 

10 79+7^** 

15 44+14.25 

< 15.84 



J033210.92- 
J033211.70- 
J033212.30- 
J033212.39- 
J033212.51- 
J033213.76- 
J033214.48- 
J033215.36- 
J033217.36- 
J033217.75- 
J033219.32- 
J033219.96- 
J033222.13- 
J033223.06- 
J033223.40- 
J033224.60- 
J033225.26- 
J033225.46- 
J033225.77- 
J033226.21- 
J033227.36- 
J033227.93- 
J033227.93- 
J033229.32- 
J033229.64- 
J033229.71- 
J033230.07- 
J033230.57- 
J033231.58- 
J033232.13- 
J033232.32- 
J033232.58- 
J033233.00- 
J033233.82- 
J033233.90- 
J033234.04- 
J033234.88- 
J033234.91- 
J033236.37- 
J033236.52- 
J033236.72- 
J033236.74- 
J033237.26- 



274722.8 
■274507.6 
274513.1 
274353.6 
274454.8 
274616.6 
274320.1 
274506.9 
274307.3 
274547.7 
274514.0 
274449.8 
274344.5 
274226.3 
274316.6 
274428.1 
274524.0 
275154.6 
274459.3 
274426.3 
275015.9 
274353.6 
275235.6 
275155.4 
274242.6 
274507.2 
274534.2 
■274518.2 
274612.7 
275105.5 
■274343.6 
■275053.9 
275030.2 
274410.0 
274237.9 
275009.7 
274440.6 
274501.9 
274543.3 
275006.4 
274406.4 
275206.9 
274610.3 



+0.39 
-0.41 
+ 1.13 
-0.08 
1 7Q+0.40 
^•'^-0.42 



1.85 
0.08 



1.35 
2.02 



1.15 
-1.35 
+0.20 

0.20 



2.86!0:« 



1 53+1-2'^ 

—^■^-1.53 



0.00+0-5^ 

2 7Q+0.76 



0.00+0-2^ 

"■•'"^-0.12 
1 oi +0.58 



2 05+° *'3 

"^■^-'-0.69 
1 O4+0.23 
^■"^^-0.23 
1 24+0.24 
^ -^^-0.25 

3.10!«J? 



62+0-25 
65 +''•5° 



i-70!^;l? 

2.9o:j:| 



o 2 S +0.53 
^•^-^-0.57 



1.27 + 0.32 

<0.31 
1.29 + 0.31 
2.44 ± 0.30 
<0.05 

< 1.97 

< 2.603 

< 1.077 
<0.28 

< 1.70 

< 1.46 
<2.17 

< 3.22 

< 1.91 
2.04 ± 0.30 

< 1.86 

< 1.49 
1.71+0.29 
0.71 ± 0.32 

<2.48 
0.79 + 0.31 

<2.02 
1.16 + 0.32 

< 1.07 
1.96 ±0.32 

< 1.61 

< 1.17 
1.67 + 0.31 

< 1.39 
2.55 ±0.31 

< 1.73 

< 1.53 
2.41+0.32 

< 2.19 
1.13 ±0.32 

< 1.64 
2.95 + 0.28 

< 2.01 
5.27 ± 0.28 
1.25 ±0.32 
2.08 + 0.23 
1.66 + 0.31 

< 1.74 



2.17' 

9 OO+0.78 
•■^°-0.92 

1.91' 



1 02+°'''* 

^•"^-0.64 

2 TT+0.53 
^■-'■^-0.57 

1.49' 
1.71+0.29 
0.71 ± 0.32 

9 79+1.13 
^•'"^-1.52 

0.79 + 0.31 

1 01 +"-2** 
0.74 + 0.22 

1 iq+0.44 
^•^^-0.47 

1.96 + 0.32 

1 OQ+0.47 
l-'^-^-0.50 

1.21+0.27 
1.45 + 0.28 
2.25 ±0.31 

0-18!if8 
1.73' 
1.53' 
2.41+0.32 
2.19' 

1 14+0.40 
^■^^-0.42 

2.95 + 0.28 

2.01' 
5.27 ± 0.28 

^•^°-0.38 

2 4Q+0.78 
"^■^^-0.94 

1.66 + 0.31 

2 4Q+0.42 
^•^^-0.44 



15.14 + 7.50 

< 16.51 
13.43+6.66 
14.95 ±4.51 

< 20.66 

< 16.66 
<9.05 

< 32.99 

< 14.56 

< 20.87 

< 20.11 

< 25.17 

< 8.83 

< 20.84 
37.95 ± 18.80 

< 8.70 

< 15.79 
32.10+15.90 
24.86 ± 12.32 

<6.92 

17.49 + 8.67 
<5.61 

8.46 + 4.19 
<7.52 
101.08 ±50.08 

< 21.08 

< 14.59 
43.69 + 21.65 

< 15.00 
20.13 ±9.98 

< 8.52 

< 16.05 
71.34 + 35.35 

< 15.79 
14.25 ± 7.06 

< 18.39 
27.90+ 13.82 

< 15.74 
12.76 ±6.32 

16.50 ±8.17 
19.17 + 9.50 

22.68 + 11.24 

< 21.02 



31.96 
24.92 



-13.75 
-13.75 
-12.07 
-12.07 

9 36+28.38 

17.49 + 8.15 

i-62!?:g 
5-05!!:^^ 

8.68!^:^^ 

J QQ 7 J +48.76 

97 74+21.80 
^'■'^-23.67 
1 S 97+6.04 
^-'■^'-6.04 

33 30+'3-^'' 

A-i 41 +20.22 
'+j.ti_2Q22 

4 64+^-^3 

< 8.48 

< 16.00 

70 0^+34.34 
'^•^■^-34.34 

< 15.80 

10 47+4-32 
^"•^'-4.51 

Q qT+6.34 
^•^-"-6.75 

28.00!||| 

< 15.76 
12.69 ± 5.23 

91 QQ + 12.66 

31 69+5 '2^ 

-^^•"^-31.69 

99 C7+IO.51 
z,z..J/_jQ5j 

CO 2'^+36.15 
-'-'•^-'-38.38 



4, 


.60 


2, 


.04 


3, 


.37 


1, 


.03 


2, 


.90 


1, 


.45 


0, 


.87 


4, 


.38 


5, 


.83 


3, 


.21 


2, 


.81 


2, 


.30 


0, 


.56 


2, 


.92 


4, 


.25 


1, 


.65 


5, 


.07 


5, 


.92 


2, 


.96 


0, 


.51 


4, 


.63 


0, 


.32 


4, 


.31 


2, 


.18 


7, 


.04 


2, 


.75 


3, 


.24 


11.19 


1 


.96 


1 


.40 





.60 


3, 


.15 


2, 


.38 


1, 


.42 


5, 


.71 


2, 


.87 





.34 


1 


.88 


0, 


,62 


4, 


,78 


1 


.17 


3, 


.05 


2, 


.44 
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Table 7. The extinction Av(Balmer), Av(lR) and the adopted Ay. The Av(IR) calculated with the lower limit of IR detection are upper limit. The 
final extinction estimated only with the limit Ljr are noted with '. Three SFR are given : SFRi^, SFRjit,, 5Fi?uv- The SFR^n, is corrected from 
extinction and aperture. To determine the unobscured 5F/fuv, we used the rest-frame 2800A luminosity from SED fitting 



Name 



Ay 

B aimer 



IR 



SFRjR 



S FRuh 



S FRuY 
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Table 8. Some important emission line ratios, the oxygen abundance in ISM and stellar mass. The uncertainties of the line ratio and metallicity are 
from uncertainties of extinction and emission line flux measurement. Galaxies in which extinction have been estimated only with AyilR^im) have 
lower limit of metallicity. When log/f23 > 1 the oxygen abundance is quoted with 
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Table 8. Some important emission line ratios, the oxygen abundance in ISM and stellar mass. The uncertainties of the line ratio and metallicity are 
from uncertainties of extinction and emission line flux measurement. Galaxies in which extinction have been estimated only with AyilR^im) have 
lower limit of metallicity. When log/f23 > 1 the oxygen abundance is quoted with 
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